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Abstract — In this paper, we characterize the throughput of 
a broadcast network with n receivers using rateless codes 
with block size K. We assume that the underlying channel 
is a Markov modulated erasure channel that is Ltd. across 
users, but can be correlated in time. We characterize the 
system throughput asymptotically in n. Specifically, we 
explicitly show how the throughput behaves for different 
values of the coding block size A' as a function of n, as 
n — > oo. For finite values of K and n, under the more 
restrictive assumption of Gilbert-Elliott channels, we are 
able to provide a lower bound on the maximum achievable 
throughput. Using simulations we show the tightness of the 
bound with respect to system parameters n and K, and 
find that its performance is significantly better than the 
previously known lower bounds. 

I. Introduction 

In this wor]^| we study the throughput of a wireless 
broadcast network with n receivers using rateless codes. 
In this broadcast network, channels between the trans- 
mitter and the receivers are modeled as packet erasure 
channels where transmitted packets may either be erased 
or successfully received. This model describes a situation 
where packets may get lost or are not decodable at the 
receiver due to a variety of factors such as channel fading, 
interference or checksum errors. We assume that the 
underlying channel is a Markov modulated packet erasure 
channel that is Ltd. across users, but can be correlated in 
time. We let 7 denote the steady state probability that a 
packet is transmitted successfully on the erasure channel. 

Instead of transmitting the broadcast data packet one 
after another through feedback and retransmissions, we 
investigate a class of coding schemes called rateless codes 
(or fountain codes). In this coding scheme, K broadcast 
packets are encoded together prior to transmission. K is 
called the coding block size. A rateless encoder views 
these K packets as K input symbols and can gener- 
ate an arbitrary number of output symbols (which we 
call coded packets) as needed until the coding block 
is decoded. Although some coded packets may get lost 
during the transmission, rateless decoder can guarantee 
that any K(l +s) coded packets can recover the original 
K packets with high probability, where e is a positive 
number that can be made arbitrarily small at the cost 
of coding complexity. Examples of rateless erasure codes 
include Raptor codes [2], LT Codes [3] and random 
linear network codes [4], where the former two are more 
efficient when K is very large and random linear network 
code is more efficient when K is relatively small and 



the field size of packets is large. The best encoding and 
decoding complexity of rateless codes (e.g. Raptor codes) 
increase linearly as the coding block size K increases. 
Further, increasing the coding block size can result in 
large delays and large receiver buffer size. Therefore, real 
systems always have an upper bound on the value of K. 

We consider broadcast traffic and a discrete time que- 
ueing model, where the numbers of packet arrivals over 
different time slots are independent and identically dis- 
tributed and the packet length is a fixed value. We let 
A denote the packet arrival rate and assume that the 
encoder waits until there are at least K packets in the 
queue and then encodes the first K of them as a single 
coding block. In this case, the largest arrival rate that 
can be stabilized is equal to the average number of 
packets that can be transmitted per slot, which we call 
the throughput. Therefore, we only need to characterize 
the throughput that can be achieved using rateless codes 
under parameters K and n. As described in Figure [T] 
the channel dynamics for the i th receiver is denoted by 
a stochastic process {Xij}j e ^, where j is the index of 
the time slot in which one packet can be transmitted 
and Xij is the channel state of i th receiver during the 
transmission of the j th packet. We capture a fairly general 
correlation structure by letting the current channel state 
be impacted by the channel states in previous I time slots, 
where I can be any number. As the number of receivers 
n approaches infinity, we show that the throughput is 
nonzero only if the coding block size K increases at least 
as fast as logn. In other words, if c = lim„_ i . 00 jj^, the 
asymptotkr] throughput is positive whenever c > 0. In 
Theorem [lj by utilizing large deviation techniques, we 
give an explicit expression for the asymptotic throughput, 
which is a function of K, n, 7 and the channel correlation 
structure. 

To study the non-asymptotic behavior of the system, 
we make a more restrictive channel assumption that the 
current channel state is impacted by only the channel 
state in previous 1 time slot, which is the so called Gilbert- 
Elliott channel model. In this case, for any finite K and 
n, we find a lower bound on the throughput in terms 
of the transmission time of a system with larger K and 
n. As a special case when the channels are memoryless, 
if is kept constant, this lower bound reveals that 
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Fig. 1: Broadcast with discrete time queueing model 



the throughput will follow a decreasing pattern as the 
number of receivers n increases. By combining this result 
with the characterization of the asymptotic throughput, 
we are able to provide a lower bound on the maximum 
achievable throughput for any finite values of K and n. 
This lower bound captures the asymptotic throughput in 
the sense that when n approaches infinity, it coincides 
with the asymptotic throughput. 

A. Related Work 

Among the works that investigate the throughput over 
erasure channels, [5]|, |]6l, J7J and ]8j are the most 
relevant to this work. In Q, the authors investigate the 
asymptotic throughput as a function of n and K and 
also show that the asymptotic throughput will be non- 
zero only if K at least scales with logn. However, they 
only consider the channel correlation model with I = 1 
and use a completely different proof technique. Moreover, 
no explicit expression on the asymptotic throughput is 
provided. In J5j and Q, two lower bounds on the 
maximum achievable rate A are provided. However, their 
bound does not converge to the asymptotic throughput 
when n approaches infinity. Moreover, our bound is shown 
to be better in a variety of simulation settings with finite 
K and n, as will be showed in Section [v| In |8|, the 
authors consider the case when instantaneous feedback is 
provided from every user after the transmission of each 
decoded packets, while we only assume that feedback is 
provided after the entire coding block has been decoded. 

B. Key Contributions 

The main contributions of this work are summarized as 
follows: 

• We give an explicit expression for the asymptotic 
throughput of the system when the number of re- 
ceivers n approaches infinity for any values of K as 
a function of n under the erasure channel with any 
levels of correlation. (Theorem [lj 



* Under the Gilbert-Elliott channel model (I = 1), for 
any finite K and n, we find a lower bound on the 
throughput in terms of the transmission time of a 
system with larger K and n. As a special case, when 
channels are memoryless (I = 0), this lower bound 
reveals that when K grows with n in a way that the 
ratio is kept constant, the throughput follows a 
decreasing pattern as n increases. (Theorem [2J 

• We provide an asymptotically tight lower bound on 
the maximum achievable throughput for any values 
of K and n under the Gilbert-Elliott channel model 
(Z = 1) and show that its performance is significantly 
better than the previously known bounds in [5] and 
@. (Theorem ^ 

The rest of this paper is organized as follows. In 
Section Jn] we describe our model and assumptions. In 
Section [Hi] we give the characterization of the asymptotic 
throughput. In Section [IV] we provide a lower bound on 
the maximum achievable throughput for any finite values 
of K and n. In Section |v| we use simulations to verify our 
theoretical results. Detailed proofs on all the theorems can 



be found in Section[VT] Finally, in Section VII we conclude 
the paper. 

II. System Model 

We consider a broadcast channel with n receivers. Time 
is slotted, and the numbers of broadcast packet arrivals 
over different time slots are i.i.d. with finite variance. We 
denote the expected number of packet arrivals per slot 
as the packet arrival rate A. The transmission starts when 
there are more than K packets waiting in the incoming 
queue intended for all the receivers. Instead of transmit- 
ting these packets one after another using feedback and 
retransmissions, we view each data packet as a symbol 
and encode the first K of them into an arbitrary number 
of coded symbols as needed using rateless code (For 
example, Raptor Code |2| or random linear network code 
[4 1) until the coding block is decoded. These K packets 
together form a single coding block with K being called 
block size. During the transmission, the coded symbols 
are transmitted one after another. 

Each receiver sends an ACK feedback signal after it has 
successfully decoded the K packets. In the following con- 
text, the term packet and symbol are used interchangeably. 

We model the broadcast channel as a slotted broadcast 
packet erasure channel where one packet can be trans- 
mitted per slot. The channel dynamics can be represented 
by a stochastic process {Xij}i<i< n j e fn, where Xij is the 
state of channel between transmitter and the i th receiver 
during the transmission of j th packet (we also call it the 
j th time slot in the i th channel), which is given by 



X^ - 



1 j th packet in the i th channel is 

successfully received 
otherwise 



We assume that the dynamics of the channels for different 
receivers are independent and identical. More precisely, 
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for all 1 < i < n, {Xy }_,->! are independent and identical 
processes. 

Since, in practice, the channel dynamics are often tem- 
porarily correlated, we investigate the situation where the 
current channel state distribution depends on the channel 
states in the preceding I time slots. More specifically, for 

T im = {Xijjj^n and fixed I, we define Him = {X im , 
Xi( m -i+i)} for m > I > 1 with Him = {0, ^} for I = 0, 
and assume that PLY l(m+1) = l|J" im ] = F[X l{m+1} = 
for al l m > Z. To put it another way, when I > 1, 



the state (Xr. 



i Xi(m-l+l) 



> I forms a Markov 



chain. Denote by IT the transition matrix of the Markov 
chain { (X im , ■ ■ • , X l{m _ l+1) ) } m>l , where 



n = [n(s, u) 



s,ue{o,i}' 5 



with 7r(s, u) being the one-step transition probability from 
state s to state u. Throughout this paper, we assume that 
IT is irreducible and aperiodic, which ensures that this 
Markov chain is ergodic [9]. Therefore, for any initial 
value Hi, the parameter ji is well defined and given by 

7l = lim F[X im = 1], 
and, from the ergodic theorem [9] we know 



lim 

m— f oo 



7=1 A y 



7i 



Since {Xij}j>i for all 1 < i < n are i.i.d., we denote 
7 = ji, for all 1 < i < n. 

Using near optimal rateless codes, such as Raptor Codes 
113, LT Codes [3] and random linear network codes Q, 
only slightly more than K coded symbols are needed to 
decode the whole coding block. For simplicity, here we 
assume that any combination of K coded symbols can 
lead to a successful decoding of the K packets. 

According to the above system model, we have the 
following definitions: 

Definition 1: The number of time slots (number of 
transmitted coded symbols) needed for user i to success- 
fully decode K packets is defined as 



The initial state for all the receivers is then denoted as 
£ = [Si, £2, ■ ■ ■ , £ n ]- 

Definition 4 (Throughput): For a system with an infi- 
nite backlog of packets, we define throughput rj(n,K) 
as the long term average number of packets that can be 
transmitted per slot. More precisely, 

ri(n,K) =Kx lim ffl, 

t-s-oo t 

where R(t) is the number of successfully transmitted 
coding blocks in t time slots. For any finite values of 
K and n, it is easy to check that {£ h ,T h (n, K)} h is a 
finite-state ergodic Markov renewal process, where £ h 
and T h (n, K) denote the initial state and the transmission 
time of the h th coding block, respectively. Then, R(t) is the 
total number of state transitions that occur in t time slots 
of this Markov renewal process. Therefore, from 11 1 Oil we 
know that 



ry(n, K) =' 



K 



K 



(1) 



E[T(n,K)} E[E[T(n,K)\£}}'' 

where the outer expectation in the last term denote the 
expectation with respect to the steady state distribution 
of the embedded Markov chain {£ h } h . 

III. Asymptotic Throughput 

Before presenting the main results, we need to intro- 
duce some necessary definitions. First, define a mapping / 
from the state space of the Markov chain {0, 1}' to {0, 1} 
as 

f((Xi m , . . . ,X( m -; + l))) = Xim- 

Then, given a real number 9, we define a matrix Ug as 

' Hs,uyf^] suem}l when/>l 
he 8 ] when / = 



Last, define a standard large deviation rate function 

A(/3,n) as 



A(/3,n) = su P {0/3-logp(n e )}, 



(2) 



Ti(K) 



j = l 



> K 



Definition 2: The number of time slots (number of 
transmitted coded symbols) needed to complete the trans- 
mission of a single coding block to all the receivers is 
defined as 



T(n,K) = max{Ti(K),i = 1,2, 



Definition 3 (Initial State): Since the current channel 
state depends on the channel states in the previous I time- 
slots, for each receiver i, by assuming that the system 
starts at time slot 1, we define the initial state of receiver 
i as 

Si = [Xi(-i + i),x i (_ l+2 ), . . . ,x i0 ] e {o,i}'. 



where p(Tle) denotes the Perron- Frobenious eigenvalue 
of Ug (See Theorem 3.1.1 in HUB ), which is the largest 
eigenvalue of II e . 

The asymptotic throughput for any values of if as a 
function of n is characterized by the theorem below: 

Theorem 1: Assume that if is a function of n and 
the value of lim 



„ ioo exists, which we denote by 



c = lim^oo j^j, then we have 



lim r](n, K) = sup < (3 



c > 



,0</3< 7 



(3) 



A(/3,n) 

Proof: see Section VI -A| ■ 
From Theorem [TJ we know that, if the coding block 
size K is set to be a function of the network size n, 
then we can characterize the asymptotic throughput when 
n approaches infinity in an explicit form. Equation ^ 
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implies that the asymptotic throughput is a function of 
7, lim n ->oo Kj log n and the channel correlation structure 
indicated by II. 

By Theorem [T] the asymptotic throughput in the special 
cases when K e o(logn) and K e w(logn) are given in 
the following corollary. 

Corollary 1.1: Assume that K is a function of n. We 
then have 

1) if K e o(logn), therQ 

lim 7y(n, K) = 0. 

n— >oo 

2) if AT e w(logn), then 

lim 77(72, X) = 7. 



Proof: 1) If if e o(logn), then c = lim r , 



and we have 



K 
log n 



C > 



-,0</3< 7 



A(/3,n) 

According to Theorem [TJ we get 

lim r)(n, K) — sup{0} 



{0}. 



2) If K e w(logn), then c 
have 



lim r , 



= 0. 

A" 



00 and we 



According to Theorem [TJ we get 

lim r)(n,K) = sup[0,7) = 7. 

n— >oo 

■ 

Corollary says that the throughput will vanish to 
as n becomes large, when K does not scale as fast as 
log n. Whereas when K scales faster than log n (Or more 
specifically when K e w(logn)), throughput approaches 
the capacity 7 of the system in the limit. It should be 
noted that Theorem [TJ together with Corollary are 
a generalized version of Theorem 1 in 0, which only 
consider the case when I = 1 and does not give an explicit 
expression for the asymptotic throughput. 

As a special case when the channels are memoryless 
(7 = 0), we can express A(/3, II) in a closed form, as shown 
in the corollary below. 

Corollary 1.2: Assume that K is function of n and the 
channels are memoryless (7 = 0), we have 

c, where c is a positive constant, then 



if limr, 



K 



lim r](n, K) - 

n— >oo 

P 



sup< ft 



log : 



7 



01 1 
- log 



^>-,0<^< 7 !>. (4) 

1 — 7 c 



Proof: When I = 0, Hg = [ye 6 ] is a degenerate matrix 
with a single entry and p(H$) — p(ye e ) = je e . Therefore 
we have, according to Equation ([2]), 

AGs, n) = log 2 + (1-0) log ^—^ 
7 1-7 



IV. Throughput lower bound for finite K and n 

For all rateless coding schemes, the encoding and de- 
coding complexity increases linearly in K, the size of 
the coding block. Moreover, the value of K determines 
the receiver buffer size. Therefore, in reality, the value 
of K is often limited by the decoder buffer size or the 
computational power of both sender and receiver. We then 
have to consider the case when K is finite and need to 
answer the following questions: For a given number of 
receivers n, channel statistics, and a maximum available 
coding block size K, what is maximum packet arrival rate 
A that can be supported by this system? For a specific 
number of receivers and channel statistics, if we are given 
a target packet arrival rate A, how can we design the value 
of K in the system such that the target arrival rate can 
be supported? 

In order to answer these questions, we make a more 
restrictive channel assumption that the current channel 
state is impacted by only the channel state in the previous 
1 time slot. In other words, we have the Gilbert-Elliott 
channel model. Under this model, for any receiver 1 < i < 
n, the channel states {AT^ j-jgN evolve with j according to 
a two state Markov chain as illustrated in Figure [2] Here, 
p i and pi are the transition probabilities between state 
and state 1, and system capacity 7 = poi/(pio +Poi)- 
Based on this model, in the theorem below, we find a 
lower bound on the throughput for any n and K in terms 
of the expected transmission time of a system with larger 
n and K. 



I-P01 P " 



1-Pu 




Fig. 2: Gilbert- Elliott Channel Model 



Theorem 2: Under the Gilbert-Elliott channel model 
(/ = 1) as illustrated in Figure [2j for any n e N K e N 
and qgN, we havq^J 

r)(n,K) > 



E[T(n a ,a(K + K ))\£ = l n „ 



where 

Kq = min < m > 



X^ 1 ~ Pw) d Pw + P01 > 1 



d=0 



Remark 2.1: Observe that K is independent of the 
choice of n and K and is only a function of channel 
dynamics. K = if and only if p i + P10 > !■ 



Proof: see Section VI-B 



When the channels are memoryless, the above theorem 
reduces to a simpler form, as shown in the following 
corollary. 



3. We use standard notations: /(n) 

if lirrin—^oo - j n ? 

9(1) 



: and /(n) 



o(g(n)) 
u(g(n)) if lim„ 



i^l diverges 

g(n) o 



4. We denote an all-one vector with dimension Til 3S 1 fj 
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Corollary 2. 1 : When the channels are memoryless (7 = 

0), for any n e N, K <E N and a e N, we have 

77(71, if) > ■q(n a ,aK). 

Remark 2.2: While Theorem Q] tells us that in order to 
achieve nonzero throughput, we can double the coding 
block size K for every quadratic increase of n, which 
is to make K/logn a fixed value, it does not tell us 
anything about how the throughput will converge as n 
approaches infinity. This corollary indicates that under the 
memoryless channel assumption, if we adapt the coding 
block size K with the increase of network size n in a way 
that K/ log n is kept as a fixed value, then the throughput 
will follow a decreasing pattern before it reaches the 
asymptotic throughput. 

Proof: The memoryless channels can be considered as 
a special case of the Gilbert-Elliott channels when p w = 
1 — 7> Poi — 7 an d K = 0, then by Theorem [2] 

^ K)> E[T{n^aK)\£ = l n ,] 

The last equation follows by noting that T(n a ,aK) is 
independent of initial state £ when the channels are 
memoryless. ■ 
By the help of the Theorem [T] and Theorem [2j we can 
get a lower bound on the maximum stable throughput 
that can be achieved for any finite values of coding block 
size K and network size n, as shown in the theorem 
below. 

Theorem 3: For a broadcast network with n receivers, 
and coding block size K, under the Gilbert-Elliott channel 
model (1 = 1), the throughput is lower bounded by 



r)(n,K) > 



K 



-n 



K + Kc 



K + Kq \ log n 
and thus the system with packet arrival rate A is stable if 

K „ {K + Kq 



A < 



K + K a 



-n 



logn 



where 



Tl (r) = sup { j3 



K 



r > 



P 



< P < 7 



m > 



E 

d=Q 



A(/?,n) 

Pio) d Pio 



and 



Poi > 1 



Remark 3. 1: As a special case, when the channels are 
memoryless, we have p w — 1 — 7 and poi = 7- It is eas Y 
to obtain that K = 0. Therefore, 



r)(n,K) 



> sup 



log 



1 



p, 1 

— log- 



P log n 

> — . < p < 7 



7 p "1—7 K 
and the system with packet arrival rate A is stable if 

P , 1 



A < sup < p 



log- 



7 



P 



P, 1 
-log 



1 



P ^ log n 

> — , < p < 7 



' 7 



A" 



Proof: From Equation ([9]) in Lemma [T] we can see 
that when AT and n are finite, the transmission time of 
a coding block T(n, AT) is light-tail distributed, meaning 
that it has finite variance. Then according to 11121 . using 
Lyapunov method we know that the queue will be stable 
if the traffic intensity of this queue, which is defined as 
the packet arrival rate A over the service rate, is less than 
1. Therefore, the queue will be stable if the arrival rate A 
satisfies 



A < sup I /i 

= v(n,K) 



K/E[T{n,K)} 



< 1 



(5) 



By Theorem [2] we know that, for any integer values of a 

aK 



r]{n,K) > 



implying that 

f](n,K) 

> lim 



A 



E[T(n a ,a(K + K ))\£ 



a{K + Ao 



c^oo K + K a E[T(n a ,a(K + K ))\£ = l n „]" (6) 

Since a(K + K ) / log n a = (K + K )/ log n for any value 
of a, then by Equation ( |17| ) in the proof of Theorem [T] 
we get 

a(K + K ) 



lim 



sup i P 
K 



E [T(n a ,a(K 
K + K Q 



\ogn 



> 



KoW 

P 



A(/3,n) 



,0</?< 7 



logn 



which, by combining Equation ^ and Equation ([6]), 
completes the proof. ■ 

In order to compare this lower bound on the maximum 
achievable rate with the existing bounds given in [5] and 
Il6l . we restate Theorem 2 in and Theorem 7 in ll6ll 
as the following. 

Theorem 4 (Theorem 2 in $5$ and Theorem 7 in /J6fj: 
In a broadcast network with n receivers, coding block 
length K and packet arrive rate A, 

1) when the channels are memoryless (7 = 0) with 
erasure probability 1 — 7 and K > 16, the system is stable 
if 

A< (l~7)if 

i£" + (logn + 0.78)vT+2.6l' 

2) For Gilbert-Elliott channels (7 = 1) with state transi- 
tion probability p w and p 01 , when 1 — p w — p 01 > and 
K > 21 log n — 4, the system is stable if 

P01K 

K + 2 v /(0.78AT + 3.37) logn + 2.61 ' 

For ease of notation let us denote the bounds given 
in Theorem [4] as the CSE bound 1 and CSE bound 2 
respectively using the initials of the authors' last name. 

Firstly we should note that the CSE bounds are only 
valid when AT and n satisfy certain conditions, while our 



A < 



6 



bound is valid for any finite values of K and n. Secondly, 
our bound converges to the asymptotic throughput in the 
sense that as n approaches infinity while keeping Kj log n 
as a constant c, our bound on the maximum achievable 
rate will converge to the asymptotic throughput with 
parameter c. Or more specifically, 

K _ ( K + Kq 



,K/5 



n=2 -p 01 =0.4,P 10 =0.4 



lim 



K\ lim 



-TZ 

K + K 

K + K 



lop 



log n 
= K{c 



lim 77(71, K), 



(7) 



which can be seen from Theorem Q] and Theorem p3] 
However, the CSE bounds are not asymptotically tight. 
When we keep the ratio K/ log n to be a constant c, as 
n or K approaches infinity, CSE bound 1 even becomes 
trivial (approach 0), which can be seen from the equation 
below. 

(l-i)K 



lim — 



(logn + 0.78)Vi^ 

(1-7) 



2.61 



= lim 

if-i-oo l + (l/ c + 0.78/K)VK + 2M/K 

-0. (8) 

Next, in Section [v| we show that our bound outperforms 
the CSE bounds under various simulation settings. 

V. Simulation 

In this Section, we conduct simulation experiments to 
verify our main results. 

A. Example 1 
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Fig. 4: Illustration for example 2 



the ratio K/\ogn, in this case, it is a constant for all K 
and is equal to the asymptotic throughput with parameter 
15/ log 2. From Figure [3] we can see that as K approaches 
infinity, the maximum achievable rate converges to our 
lower bound (which is also the asymptotic throughput 
in this case) in a decreasing manner, which validates 
Theorem [T] Corollary 2.1 and Theorem [3] 

In this case, we also plot the CSE lower bounds given 
by Theorem 4. From the figure we can see that our 
bounds outperforms the CSE lower bounds. CSE bound 
1 gradually approaches zero as indicated by Equation Q 
while our bound is a constant value and asymptotically 
tight as shown in Equation ([7]). 
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This example verifies Theorem [T] Corollary 2.1 and 
Theorem [3] We choose a memoryless channel with 7 = 
0.5. By keeping Kj log n as a constant 15/ log 2, we change 
K from 5 to 300 and calculate the maximum achievable 
rate, which is 77(77,, K), through simulations for each pair 
of (K,n). Since the value of our bound is a function of 



B. Example 2 

In order to verify Theorem [3] under the Gilbert-Elliott 
channel model, we choose the state transition probability 
Pw = P01 = 0.4. It is easy to obtain that 7 = 0.5 and 
K = 1. By keeping K/ log n as a constant 5/ log 2 and 
changing K from 5 to 100, we plot in Figure [4] both the 
simulation result of the maximum achievable rate and our 
lower bound shown in Theorem [2] Again, as we can see 
from the figure, the lower bound becomes tighter as K 
increases. Eventually the maximum achievable rate will 
converge to the lower bound as shown in Equation ([T]). 
However, neither of the CSE lower bounds is valid under 
this system setting. 

C. Example 3 

In this example, we conduct three set of experiments 
under the memoryless channel assumption with different 
values of AT as a function of n, and show that our 
bound outperforms the CSE bound in all these simulation 
settings. 

In the first case, we set the coding block size K to be the 
same as the network size n and change n from 5 to 300. 
We plot the simulation result of the maximum achievable 
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Fig. 5: Illustration for example 3 
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rate as well as our bound and the CSE bound in Figure [5a] 
since in this case K scales faster than log n, the achievable 
rate will approach system capacity 7 as the network size 
n grows. 

In the second case, we assume that the number of 
receivers is fixed to be 10 and we increase coding block 
size K from 5 to 300. The simulations result, together with 



the two bounds, are plotted in Figure 5b In this case, the 
achievable rate will also approach system capacity 7 as n 
increases. 

In the final case, as shown in Figure [5cJ we keep the 
coding block size to be a constant 80 and increase the 
number of receivers from 5 to 100. Since K does not 
increase with logn at all, the achievable rate will vanish 
to as n grows. 

From Figures [3] [5aJ [5b] and [5c] we can see that our 
bound obtained by Theorem [3] is significantly better than 
the lower bounds achieved in |]5]| and |[6]| in all these four 
different cases. 



VI. Proofs 

A. Proof of Theorem [7] 

In order to prove Theorem [T] we first need the following 
lemmas (Lemma [TJ Lemma [2] and Lemma [3]). 

Lemma 1: For any /3 e (0, 1) and any values of £ , we 
have 



T{n,K) > 



K = k,£ 



-l- (l- e -|A(/3Wl(/3<7)+ff(/3,fc,£)y l 



(9) 



where 



g(/3,k,£) e 



o{k) as k — > 00 if /3 < 7 
o(l) as k —> 00 if > 7 



Proof of Lemma [I] From definition ([!]) and ([2]), we 

have, for any i, 

n 

{T{n, K) < t, £ } = p| {T^K) < t, £i} . 

i=l 



Therefore, we have 



¥[T(n,K) > t\K = k,£] 
=1 -¥[T(n,K) < t\K = k,£] 

n 

=l-l[(l-P[T i (K)>t\K = k,£ i \) 



(10) 



Let t = 4, from definition [l] we can get, for any 1 < i < n, 



Tm > - p 



K = k,£i 



2^=1 A ij 



k/0 



and 



lim 

k— >oo 



logP Ti(K) > 1 



K = k,£i 



k/P 



= -A(/3,n)l(/3< 7 ), 
(11) 



with the last equation being a direct application of 
Theorem 3.1.2 in HI 1 B (Gartner-Ellis Theorem for finite 
state Markov chains). Notice that the right hand side of 



Equation (111 is fixed for all possible values of i and 
as long as the values of {3 and II are fixed. Then the proof 



completes by combining (101 and (111. 



Lemma 2: Assume A; is a function of n and denote k := 
k(n), and define f(k,(3,£) := e iHP,n W <i)-9(PX£) > 
then we have 



1) For a fixed (3 e (0, 1), if lim^oo /(fc( „" 

k(n) 



0, then 



lim I 

n— f 00 



T(n, K) > 



K = k{n),£ 



= 0. (12) 



2) For a fixed (3 e (0, 1), if lim r , 
then 



f(k(n),f3,£) 



CO, 



lim I 

n— >-oo 



T(n, K) > 



k(n) 



K = k(ri),£ 



= 1. (13) 
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Proof of Lemma |2j According to Lemma [T] and the 
definition of f(k(n), (3,£), we have 



T(n, K) > 



k(n) 



K = k{n),£ 



f(k(n),(3,£) 



f(k(n),f3,£) 



Since the function (l — i) with domain (l,+oo) is 
a bounded and strictly increasing function with region 

(0, e" 1 ) and the fact that f{k,f3) > 1, we know that if 
lim„ 



lim inf I 

n— >oo 



T{n,K) > 



oo, then 

k(n 



=1 — lim sup 



1 



1 



K = k(ri),£ 

/(fc(n),/9,£)' 



f(k(n),(3,£) 



/(*(»), 



> 1 - lim sup e fW").f>,e) 
=1, 



which, together with the fact that P[T(n, X) > ^ |if = 



k(n),£] < 1, yields Equation (131. 

If lim„^no ttutz\ a c\ = 0, then f{k{n),f3,£) — > oo as 



J n->oo f(k(n),p,£) ~ 

n-Kxi, which results in 



lim ( 1 

n— >oo 



1 



f(k(n),P,£ 
Then we can obtain 



f(k(n),j3,£) 



lim sup I 



T(n,#) > 



=1 — lim inf 

n— fco 

= 1- 



j8 
1 



K = k(n),£ 

f(k(n),P,£)' 



lim inf 



lim 1 - 



1 



f(k(n),(3,£) 



f(k(n),p,ey 



/(fc(n),|9 ( 6) 



=1-1 = 0, 

which leads to Equation ( |12p . ■ 
Lemma 3: Let {/i„(a:)} be a set of Lebesgue measur- 
able functions defined on [0, oo) and /i„(x) converges to 
l(x < y) almost everywhere for some y > 0. If /i„(a;) is 
a decreasing function of x and have the range [0, 1] for 
any n e N, then h n {x) converges globally in measure to 
l(x < y). 

Proof: Choose e > 0. Since h n (x) converges to l(x < 
y) almost everywhere, for any 5 > 0, we can find N £ N 
such that for any n > N, we have 

\hn (y-S/2)-l\<e 
\h n (y + S/2)-0\<e. 



Since < h n (x) < 1 for any x e [0, oo) and h n {x) is a 
decreasing function of x, we know that, for any n > M, 

h n (x)>l-e Vx<y-5/2 
h n (x) < e Vx > y + <V 2 - 

Therefore, for any n > N, 

u({\h n [x) - l{x < y)\ > e}) 
< 2 y([y-5/2, 2/ ])+K[y,y + V2])=5, 

where z/ is the Lebesgue measure. Since e and 8 are 
arbitrarily chosen, from the above inequality we know 
that h n (x) converges globally in measure to l(x < y). 

■ 

With Lemma [Tj2land|3] established, we now turn to the 
proof of Theorem [T] 

Proof of Theorem [I] Since K is assumed to be a 
function of n, we denote this function as k(n). According 
to definition [4] we have, 



lim (j](n,K)) 

n— >oo 

= lim E 



E [T(n,K)\£] 
K 



(14) 



Next, we obtain the value of lim^oo E[T(n, K)\£]/K and 
show that it is independent of £. Note that 

~T(n, K) 



lim E 

n— ¥oo 



= lim 

n— >oc 

— lim 



K 



£ 



[T(n, K) > s\K = k(n),£] 
k(n) 



ds 



' [T(n, K) > k(n)u\K = k(n),£] du. (15) 



According to the assumption that lim n ^oo k(n)/\og(n) = 
c, we have 

lim 

n— f oo 



e e 



l -A( / 3,n)l(/3< 7 )- s (/3,fc,£) 



U c> 

oo c < 





A(0,n)l(/3< 7 ) 

g 

A(^,n)i(^< 7 ) 



Since 
and 



A(,3,n)l(/3<7) I/ 3 



=o = 0, lim 



/3^7 A(g,njl(/3<7) 



Ap n)i(;3<7) ^ s a mon °tone increasing function on 
the domain (0,7), the equation c = A ^ n ) / i(^<7) nas on ^ 
one solution of /3, which we denote as 



f3 c = sup <^ f3 



c > 



,0<^< 7 



A(/?,n) 

Then, by Lemma [2j we get 

lim P [T(n, K) > k(n)u\K = k(n),£ 



1 if u < j- 
if u > £ 



(16) 



We let h n (u) = W[T(n,K) > k{n)u\K = k(n),£]. Equa- 
tion (16 1 implies that h n (u) converges to l(u < l/(3 c ) 



pointwisely. Since h n (u) is a decreasing function of u and 
has the range [0, 1] for all n, by Lemma [3] we know that 
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h n {u) globally converges in measure to 1(« < l//? c )- We 
also know that the set of function {h n (u)} is uniformly 
bounded. Then we can apply Vitali convergence theorem 



to Equation (151 to exchange the limit and integral and 
obtain 



lim E 



T{n,K) 



K 



lim ] 

n— >oc 



1 



[T(n,K) > k{n)u\K = k{n),£\ du = — . 

Pc 

(17) 



Note that the above result is independent of the choice of 
the initial state £. Since the cardinality of the state space 
of £ is finite for a finite value of n, we can exchange 



the limit and expectation in Equation (14), which, after 
combining with the above equation, completes the proof. 



B. Proof of Theorem [2] 

Let us define two random variables and under 
the Gilbert-Elliott channels as 



J"(°) = min < m 

m 

= min < m 



E*lj > 1,^10=0 



> 1,^10 = 1 



(18) 
(19) 



In order to prove Theorem [2j we first need the following 
lemma. 

Lemma 4: Let {T^} den be i.i.d. random variables 
with the same distribution as then we have 

K + l 



meaning that Y^dli is stochastically greater than or 
equal to T(°\ where 



Kq = min < m > 



E( x ~Pw) d Pw +P01 > 1 



(20) 



d=0 



Proof of Lemma [4} First observe that X)dLoQ ~ 
Pw) d Pio = 1> which makes sure that K in Equation (201 
is well defined. 

Then according to the definition of and in 



Equations (18) and (191, we have, for any integer 1 < 

t<K + l, 



K + l 



E r r>* 



1 > 



T (o) > t 



(21) 



and for any integer t > K + 1, 

'K a + 1 



E p [t {1} > 1 

d=l 
K +l 

> E p F (1) > 1 

d=l 
K +l 

= E pioCi-pio^-'a-Poi) 



p 


T (l) = { 


P 


r<°> > t - d 


p 


T (i) = { 


d-l 

P 


T (o) > t _ ! 



>(l-poi)(l-poi)' 



T (0) > t 



with the last inequality followed by the definition of 
K in Equation ( |20[ ). The above equation, together with 
Equation ( |21| ), completes the proof. ■ 
With Lemma [4] established, we now turn to the proof 
of Theorem [2] 

Proof of Theorem [2j Under the Gilbert-Elliott channel 
assumption as illustrated in Figure [2] let T(n, K, £) be 
defined as T(n, K) with initial status E. Then according 
to Definitions [T] and [2] and Equations ( |18p and ( |19P we 
can express T{n,K,E) as 

where {T^}^- s n are i.i.d. random variables with the 
same distribution as T^ a \ and {7y^}i,.jeN are i-i-d- ran- 
dom variables with the same distribution as Simi- 
larly we can express T(n a ,aK 1 1„ Q ) and T(n a , aK, 1„«) 
as 

r(n,x,l„)=max n ||j7;w|, (22) 



T{n a ,aK, l r 



max 

Ki<n° 



E^ 



(23) 



First, by Lemma [4] we know that, for any 1 < i < n and 
any initial status £ , 



K + l K+K K 

E T ff } + E ^^4 £i) +E^ 

J=l j=Ka+2 J=2 

implying that 

which yields 

E[T(n,K + K ,l n )] >E[T(n,K,£)]. (24) 
Next, we will show that 

E[T(n a ,aK,l n c,)} > aE[T(n,K,l n )]. 



10 



Let us deonte 



rK 



ST 



E 



T. 



(i) 



j = l + (r-l)if 



Then we know that {S[} iS N.reN are i-i-d- random vari- 
ables. Equation (\22ty and (\23ty can be rewritten as 



(25) 
(26) 



T(n, K, l n ) = max 5** 

l<2<n 



T(n Q ,aX,l„o) = max VS[. 



Instead of viewing Equation (261 as a 1-dimensional 
maximization over n a points, we can think of it as an a- 
dimensional maximization over n a points where we can 
choose a coordinate from 1 to n on each dimension and 



therefore can further rewrite Equation ( 26 1 as 

T(n a ,aK,l n a) = max max ... max ST. 



r=l 



i 2 ,■■■,ia) , 

(27) 



where 



s; 



(ii,i2, —,<<») 



S- 



and i„ can be viewed as the coordinate in the u th dimen- 
sion. 



Next, we use Equation (271 to build a lower bound on 



the expection of T(n a ,aK, l n «). 

For fixed values of is, . . . ,i a , let us find a i* such 
that 

il(i 2 ,...,i a ) = arg max S ( \ 4 4 ); (28) 

l<ii<n v ' 1 ' ' 

which we denote as i\ for short. Then according 
to Equation (271, we can find a lower bound for 

E[T(n a ,aK, ln a )] by choosing i x = i£ (i 2 , . . . , i a ) for all 
possible values of i^, 13, ■ ■ ■ , i a > which is 



E[T(n a ,aK),l 
=E 



max max . . . max 

l<zi<nl<22<n l<i a <n 



(a) 
>E 



=E 



max . . . max 

l<-i 2 <n l<i a <n 



max . . . max 

l<-i 2 <n l<i a <n 



r=l 

r=l 
/ a 

^ / S f (i*,i 2 ,-- ■ 
\r=2 



i 2 ,...,i a ) 



(29) 



Since the choice of ij is only sub-optimal, the inequality 
(a) in Equation <\29} should be strict inequality. No- 



tice that according to Equation (281, for any values of 

i 2 , hi - ■ ■ , ia, we have 



which, combining Equation (251 and the fact that {ST} 
are i.i.d. random variables, yields 



E 



°(il,i 2 ,...,i a ) 



=E 



max S t , 



l<ii <n 

=E 
=E [T(n, K, 1 



max S} 

Ki<n 



(30) 



As a second step, for any values of i 3 , 14, ... ,i a , let us 
, as 

il(il,i 3 ,...,i a ) = arg max Sl, i i y 



define i 2 as 



Then similarly as Equation ( [29] ), by fixing i 2 to be i 2 , we 
can obtain 

E[T(n a ,aK,l na )] 
> E max . . . max 

l<i3<n 1 < 2- a < n 



ic) 



\r=3 



Also, for any values of z 3 , 14 . . . , i a , we have 



E 



max SL ■ n 



=E 

=E [T(n,K, 1„)]. 
By defining , . . . , z* in a similar way 

*u(*1j • • • J li • ■ • J *ct) 

= arg max S^-» , ■ > 

and iterating the above step, we can get 

E[T(n a ,aK,l na )] 
> E [ S lit,H,-^) + S (iir 2 ,-,i^) + •••+# 



(31) 



r=l 



( = 5 oE [T(n, AT, 1„) 



(32) 



Equation (b) follows from the fact that 
{S(i* i* i*)}i<r<a are independent random variables 



and equation (c) follows from Equations (30), (31), and 
iterative steps. By combining Equations Q, (241, and 
(|32|), we have, 



77(71, A") = 



> 



A" 

E[E[T(n,AT,£)]] 
A" 

E[T(n,A- + X ,l„)] 
aAT 



> E[T(n°,a(K + K ),l na )} 
K a(K + A' ) 

~ K + K a E [T(n<*,a(K + K ), l n » )] : 

which completes the proof. 
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VII. Conclusion 

In this paper, we characterize the throughput of a 
broadcast network using rateless codes. The broadcast 
channels are modeled by Markov modulated packet era- 
sure channels, where the packet can either be erased or 
successfully received and for each receiver the current 
channel state distribution depends on the channel states 
in previous I packet transmissions. 

We first characterize the asymptotic throughput of the 
system when n approaches infinity for any values of 
the coding block size K as a function of number of 
receivers n in an explicit form. We show that as long as 
K scales at least as fast as logn, we can achieve a non- 
zero asymptotic throughput. Under the more restrictive 
Gilbert-Elliott channel model (7 = 1), we study the case 
when K and n are finite. For any K and n, we find a lower 
bound on the throughput in terms of the transmission 
time of a system with larger K and n. As a special case 
when channels are memoryless, this result shows that, by 
keeping the ratio if/ log rt to be a constant, the system 
throughput will converge to the asymptotic throughput 
in a decreasing manner as n grows. By the help of these 
results, under the Gilbert-Elliott channel model, we are 
able to give a lower bound on the maximum achievable 
throughput (maximum achievable rate), which is a func- 
tion of K, n and state transition probabilities p i and p w . 
In contrast to the state-of-the-art, we analytically show 
that our bound is asymptotically tight when Kj log n is 
fixed as n approaches infinity. Further, through numerical 
evaluations, we show that our bound is significantly better 
than existing results. 
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